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Abstract. We analyze a morphologically-selected com- 
plete sample of 52 late-type (spiral and irregular) galaxies 
in the Hubble Deep Field North with total K-magnitudes 
brighter than K=20.47 and typical redshifts z ~ 0.5 to 
1.4. This sample exploits in particular the ultimate imag- 
ing quality achieved by HST in this field, allowing us to 
clearly disentangle the early- from late-type galaxy mor- 
phologies, based on accurate profiles of the surface bright- 
ness distributions. Our purpose was to investigate system- 
atic differences between the two classes, as for colours, red- 
shift distributions and ages of the dominant stellar popu- 
lations. Our analysis makes also use of an exhaustive set 
of modellistic spectra accounting for a variety of physical 
and geometrical situations for the stellar populations, the 
dusty Interstellar Medium (ISM), and relative assemblies. 
The high photometric quality and wide spectral cover- 
age allow to estimate accurate photometric redshifts for 
16 objects lacking a spectroscopic measurement, and al- 
low a careful evaluation of all systematics of the selection 
[e.g. that due to the surface- brightness limit]. This sam- 
ple appears to miss significantly galaxies above z = 1.4 
(in a similar way as an early-type galaxy sample previ- 
ously studied by us), a fact which may be explained as a 
global decline of the underlying mass function for galax- 
ies at these high redshifts. Differences between early- and 
late-types are apparent - particularly in the colour dis- 
tributions and the evolutionary star-formation (SF) rates 
per unit volume -, although the complication in spectro- 
photometric modelling introduced by dust-extinction in 
the gas-rich systems prevents us to reach conclusive re- 
sults on the single sources (only future long- wavelength IR 
observations will allow to break the age/ extinct ion degen- 
eracy). However, we find that an integrated quantity like 
the comoving star-formation rate density as a function of 
redshift is much less affected by these uncertainties: 
by combining this with the previously studied early-type 
galaxy sample, we find a shallower dependence of ^(z) on 



z between z — 0.2 and z = 1.5 than found by Lilly et al. 
( |l995| ). Our present results, based on a careful modelling 
of the UV-optical-NIR SED of a complete galaxy sample 
- exploiting the observed time-dependent baryonic mass 
function in stars as a constraint and attempting a first- 
order correction for dust extinction - support a revision 
of the Lilly-Madau plot at low-redshifts for both UV- and 
K-band selected sa mples , as suggested by independent au- 
thors (Cowie et al. |1999[). 
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1. INTRODUCTION 

Cosmogonic models, in particular those based on the hi- 
erarchical clustering of cold dark matter halos, now in- 
cluding detailed physical descriptions of gas cooling, star- 
formation and feed-back processes in the baryonic com- 
ponent, make specific predictions about the evolutionary 
history of galaxy populations as a function of their mor- 
phology. 

Basically, in the hierarchical scheme, forming galaxies 
acquire angular momentum from tidal interactions with 
the surrounding structure and then dissipate and collapse 
preferentially along the rotation vector and tend to form 
flattened rotational-supported structures (disk galaxies). 
There are indeed indications of a substantial population 
of large structures of this kind up to the highest red- 
shifts from absorption-line studies in the distant quasar 
spectra (e.g. Wolfe A.M. |1999[ ). The pressure-supported 
stellar bulges dominating E/SO galaxies, in this scheme, 
would originate from the violent relaxation and dynami- 
cal evolution following strong interactions and mergers of 
primordial disk galaxies. At the zero-th order, these mod- 
els predict that spheroidal galaxies are assembled some- 
what later than spirals, although their stellar populations 
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might not differ much if the merger occurs among gas- 
poor systems. One such extreme case has been discussed 
by Kauffmann & Chariot (1998) based on the 0=1 stan- 



dard CDM cosmology, predicting a substantial dearth of 
spheroidal galaxies already by z=l, most of them being 
formed at lower z. 

More recently it has been pointed out that, within this 
scheme, the morphological appearence of a galaxy may 
repeatedly change with cosmic time not only from a late- 
to an early-type following a merging event, but also from a 
spheroidal to a disk configuration following the acquisition 
of new infalling gas from the environment (e.g. Ellis R.S. 



1997). This may reflect in long- wavelength (V to K) colour 
distributions un-distinguishable between early- and late- 
types, while shorter-wavelength (U to V) colours would be 
dominated by the on-going SF in disks. 

An opposite pattern is contemplated by the "tradi- 
tional" models of galaxy formation, assuming that massive 
galaxies, in particular elliptical and SO's, originated first 
at high redshifts as single entities from rapid homologous 
collapse of primordial gas. Gas-rich systems, in this view, 
form instead more quiescently from progressive inflow of 
gas into the dark matter halos during most of the Hub- 
ble time. This formation scenario then predicts a marked 
differentiation in colours and ages for the stellar popula- 
tions of the two classes of galaxies, late-type galaxies con- 
taining much younger stellar populations on average than 
the early-types. Also a substantial population of massive 
spheroids would be expected to be visible at z ~ 2 in this 
case. 

The ultra-deep integrations at various wavelengths 
performed by the Hubble Space Telescope in the Hubble 
Deep Field (Williams et al. 



1996, we consider here only 



the survey in the North area) offer an extremely valu- 
able dataset to study morphological properties of high- 
redshift galaxies. Furthermore, the very accurate photom- 
etry achievable in such deep images allows accurate es- 
timates of the photometric redshifts for vast numbers of 
faint galaxies in the field. 

We have recently exploited these data to study the 
colours, masses, age distributions, and the star-formation 
history of a sample of elliptical-SO galaxies (Franceschini 
et al. 1998, FA98 hereafter). The basic result was to find 



colours indicative of wide ranges of ages for the stellar pop- 
ulations and a remarkable absence of objects at z > 1.3, 
both facts telling against the predictions of the "tradi- 
tional" monolithic formation scenario. 

As a natural complement, we present in this paper an 
analysis of late- type and irregular galaxies in the HDF. 
Similarly to what we did there, our primary selection is 
in the K-band, obtained from a deep KPNO image, to 
minimize the biases in the sample due to the effects of 
K- and evolutionary corrections. The completion of our 
previous analysis of E/SO to account for the complemen- 
tary set of late-type systems is also needed for a global 
evaluation of the star-formation history as a function of 



redshift. The advantage of our approach over previous at- 
tempts (Lilly et al. 1996, Cowie et al. 1999) is in our care- 
ful treatment of dust extinction from a detailed fitting of 
the UV-optical-NIR spectral energy distributions (SED). 
In addition, the detailed knlowledge of the near-IR (NIR) 
spectrum for sources at the relevant redshifts is informa- 
tive on the baryonic mass function in stars, which provides 
an essential constraint on the cumulative star-formation 
rate as a function of time. 

In Section 2 we discuss the selection scheme and pho- 
tometric corrections used to construct a complete K-band 
flux limited sample of late-type galaxies. In Section 3 
we describe the population synthesis code that we used 
to model the optical-NIR SEDs of our sample objects, 
taking into full account the effects of a dusty interstellar 
medium in the galaxy spectra. Our main results are then 
reported in Section 4, where we perform detailed anal- 
yses of the space distributions, colours and ages of the 
stellar populations of field spiral and irregular galaxies, 
compared with ellipticals and SO. We discuss the difficul- 
ties inherent in the spectral modelling of gas-rich systems 
affected by dust extinction. We finally attempt to con- 
struct the global star formation histories of field galaxies 
(E/SO+spiral/irregulars). In Section 5 we summarize our 
main conclusions. 

We anticipate that the results of the present analysis, 
based on a survey over a very small sky field, are to be 
considered as only tentative, untill larger areas will be 
surveyed to similarly deep limits. 

We adopt Ho = 50 Km s -1 Mpc^ 1 throughout the 
paper. For consistency with FA98 the analysis is made 
assuming go = 0.5, and zero cosmological constant A. 

2. SAMPLE SELECTION AND PHOTOMETRY 

The Hubble Deep Field North has been observed in 4 
broad bands (F300W, F450W, F606W, F814W) for a to- 
tal of 150 HST orbits by Williams et al. (1996), and con- 
stitutes the deepest ever exposure on a small sky area. 
Accurate photometric data in the four bands have been 
published for thousands of faint galaxies by the authors. 

Dickinson et al. ( |1997| ) observed the HDF-North in 
the near-IR with the IRIM camera on the KPNO 4 m 
telescope. The camera employs a 256 x 256 NICMOS-3 
array with 0" .16 pixel -1 , but the released images were ge- 
ometrically transformed and rebinned into a 1024 x 1024 
format. IRIM exposures have been secured in the J, H 
and K filters, for a total of 12, 11.5 and 23 hours, respec- 
tively. Formal 5 a limiting magnitudes for the HDF/IRIM 
images, computed from the measured sky noise within a 
2" diameter circular aperture, are 23.45 mag at </, 22.29 
mag at H, and 21.92 mag at K, whereas the image quality 
is - 1".0 FWHM. 

Our sample of galaxies has been extracted from the 
HDF/IRIM if-band image through a preliminary selec- 
tion based on the automatic photometry provided by SEx- 
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tractor (Bertin & Arnouts 1996). It is flux limited in the 
K band and it excludes early-type galaxies, i.e. objects 
whose surface brightness distribution is dominated by a 
de Vaucouleurs profile. 

To determine the limit of completeness in the K band 
(hereafter Kl) for inclusion in our sample, we followed 
the same empirical procedure described in FA98. Briefly, 
a large number of toy galaxies with exponential profiles 
(produced with the IRAF-MKOBJECTS tool) have been 
used to check the performances of Sextractor in estimating 
the K-band magnitudes {Ksex) of late- type galaxies. This 
allowed us to determine the magnitude Ku m below which 
the scatter aK SEm of the measured magnitudes turns out 
to be lower than some given value a ma x (the vertical line 
in Figure la corresponds to K\ im — 21 and <r max — 0.18). 
Moreover, using only galaxies with < K$Ex >< Ku m , we 
have derived the following empirical relation between the 
bias AK =< K$ex > —K trU e and the effective surface 
brightness < pt^(SEx) > (see Figure lb): 

log AK = -1.059 + 0.188 < /if > +0.061 < /if > 2 

+0.006 < /if > 3 , (1) 

where < /if >=< fi^(SEx) > -21. This relation pro- 
vides the true total K magnitudes from the SExtrac- 
tor measured flux. Then simulations have shown that 
the K-band image has a surface brightness limit of < 
H?(SEx) >= 23. 

By analogy with FA98, we have first produced a cat- 
alog of morphologically selected late-type objects with 
Ksex < Kum = 21.0; then we have used effective ra- 
dius estimates from high resolution (HST) optical imag- 
ing to derive the effective surface brightness of galaxies 
(< /tf >); finally, we have applied to the K$Ex magni- 
tudes the statistical corrections given by equation [1] and 
we have included in the final sample only galaxies with 
corrected magnitudes less than or equal to Kl = Ku m — 
AK max - a ma x = 20.47 mag (we assume AK max ~ 0.35 
from Figure lb). 

A total of 176 objects with K$Ex < 21 mag were de- 
tected by SExtractor in the IRIM if-band image. After 
carefull inspection of the high resolution HST images, we 
rejected all elliptical and SO galaxies (including the 35 
early-type objects identified in FA98) and the stars. A few 
objects were also rejected from the sample due to their po- 
sition in the frame (at the edges of the image the noise is 
higher and the magnitude estimate is likely to be uncer- 
tain). We then produced a first preliminary, incomplete 
sample of late-type galaxies. 

The effective radii r e were estimated running SEx- 
tractor on the WFPC2 Vqqq frame with the parame- 
ter FLUX -RADIUS, providing the radius containing 
half of the total emitted flux. The surface brightness 
< fif-(SEx) > was evaluated for each galaxy and the sta- 
tistical corrections AK were computed using the eq (1). 

The final complete sample of late-type galaxies with 
Kcorr < 20.47 consists of 52 objects over the HDF area 
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Fig. 1. Left panel: standard deviation of the SExtractor 
magnitude estimates as a function of the average magni- 
tude for galaxies in the simulated images: the standard de- 
viation is very small for < Ksex >< 21™. 0. Right panel: 
difference between true flux and the SExtractor flux as a 
function of the average surface brightness for the subsam- 
ple of simulated galaxies with < Ksex >< 21 m .O. 



of 5.7 square arcmin. For 36 objects we have the spec 
troscopic redshift (Cohen et al. 1996, Cowie et al 



1999 



Fernandez- Soto et al. 1998| ), while for the remaining 16 we 
measured it from our photometric analysis as described 
below. 

A procedure analogous to that outlined for the K band 
was used to derive, for each object of the selected sample, 
the corrected magnitudes in the J and H bands. 

The optical magnitudes in the F300W, F450W, 
F606W and F814W bands (U, B, V and I in Table 1, re- 
spectively) have been computed again with SExtractor on 
the high resolution WFPC2 images (no corrections being 
applied in this case). Magnitudes are in th e AB system, 
defined by the relation (Oke & Gunn 1983): 



AB = -2.5 log F u -48.60 



(2) 



where F v is the flux in erg cm~ 2 Hz -1 , the constant bee- 
ing choosen so that AB = V for an object with flat spec- 
trum. 

Some data on the sample are listed in Table 1 . Column 
1: our identification; column 2-4: coordinates a and 5 (at 
J2000). To these must be added 12 hours 36 minutes (RA) 
and 62 degrees (Dec); column 5: r e is the effective radius, 
in arcsec, derived from HST images; column 6-9: optical U, 
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Fig. 2. Observed broadband spectra for sixteen galaxies in our sample, fitted with the synthetic models described in 
Section 3. 



B, V, I magnitudes in the AB system (see above); column 
10-12: near-infrared J, H, K corrected magnitudes in the 
standard system; column 13: redshift of each object. Val- 
ues in brackets are photometric redshifts, while the other 
are all spectroscopic. 



stellar medium in galaxy spectra. We defer the reader to 
that paper for a through description of this model and for 
precise definitions of the parameter, while for convenience 
we summarize the main features below. 



3. MODELLING GALAXY SEDs IN THE 
PRESENCE OF A DUSTY ISM 

The optical-NIR SEDs of our sample objects have been 
modelled u sing t he population synthesis code GRASIL 
(Silva et al. 1998 ), taking into full account the effects (op- 
tical extinction and thermal reprocessing) of a dusty inter- 



3.1. The GRASIL code 

The code provides a self-consistent description of the for- 
mation and evolution of a galactic system in its various 
stellar and ISM components, including its secular evolu- 
tion during the Hubble time and episods of enhanced star- 
formation possibly following interactions and mergers. 
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As a preliminary step the code allows to solve the equa- 
tions ruling the chemical evolution, providing the star for- 
mation and metallicity histories SFR(t) and Z(t) as a func- 
tion of time. The computations presented here were per- 
formed adopting one-zone (no spatial dependence) open 
models including the infall of primordial gas, according to 
the standard equations of galactic chemical evolution. As 
usual, the star formation rate is determined by the amount 
of gas in the system according to a Schmidt-type law [] 

SFR(t) = vM g [tf. 

We have generated 3 different SFR(t), in order to provide 
a wide range of spectral evolution patterns. The peak oc- 
curs at about 1, 2 and 3 Gyr (hereafter model (a), (b) 
and (c) respectively), getting broader from (a) to (c). As 
a result, half of the final stellar mass M (i.e. at 13 Gyr) 
has been assembled at galactic times of 2, 3.7 and 4.7 Gyr 
in the three cases respectively. A standard Salpeter IMF 
between 0.1 and 100 M Q is assumed. 

As described by Silva et al. (1998), GRASIL calcu- 
lates self-consistently the absorption of starlight by dust, 
the heating and thermal emission of dust grains, for an 
assumed geometrical distribution of the stars and dust, 
and a specific grain model. 

In the GRASIL model several parameters affect the 
overall modifications imprinted by dust on the SED. How- 
ever, if we confine ourselves to the attenuation of stellar ra- 
diation in the optical/UV/NIR bands, we can obtain most 
of the possible spectral behaviours by adjusting only two 
quantities: the escape timescale t (see Eq [8] in Silva et 
al. 1998 for a precise definition) of newly formed stars from 
parent molecular clouds (MCs) and the total mass of dust. 
Indeed t Q controls the fraction of light from very young 
stellar generations hidden inside MCs and converted to 
IR photons, since the MCs optical thickness is very high 
below ~ l^tm (cfr. Silva et al. 1998). On the other hand, 
the effects of the diffuse (cirrus) dust depend on several 
quantities: the radial and vertical scale lengths for stars 
and dust distributions Rd and z d , the residual gas in the 
galaxy M g , the dust to gas ratio 5 and the fraction of gas 
which is in the MCs component f mc . However we found 
that most, if not all, the possible attenuation laws of the 
diffuse dust, arising from different choices of these quan- 
tities, can be closely mimicked by simply adjusting the 
amount of gas, while fixing the other quantities to the 'typ- 
ical' values: R d = 3.5(M/10 11 M Q ) 1 / 3 Kpc, z d = 0.1 R d , 
S = 0.01 and f mc = 0.5. Obviously, while different choices 
of M g , R dl Zd, S and f mc can yield similar attenuation 
laws on the optical spectrum, the spectral shapes of the 
corresponding IR continuum re-radiation can be rather 
different. 

Strictly speaking the residual gas M g is not a param- 
eter, being instead the outcome of the chemical evolution 
code, through the Schmidt law. However we use the trick 




1 The code can be downloaded from http://grana.pd. astro. i\ 
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Fig. 3. Comparison of photometric redshifts, based on 
seven-band spectral data, with spectroscopic redshifts. Er- 
ror bars refer to interval solutions with more than 90% 
confidence level, derived from % 2 fitting using models [a] 
and fcl. 



of forcing M g to different values, in order to describe with 
a monoparametric sequence the effects of a global attenu- 
ation on the SED. Besides this, a larger 'freedom' on M g 
takes into account that the Schmidt law should not be 
taken too literally, as a strict relationship between the to- 
tal gas content and the SFR in the system. The law may 
only provide an order of magnitude description, in par- 
ticular for the secular evolution of the SFR, the so-called 
"inactive phase" of galaxy evolution bringing essentially 
to the formation of spiral disks. Several other physical 
parameters influence the rate of star-formation with re- 
spect to the simple available amount of residual gas, in 
particular the gas pressure and temperature, which may 
drastically change as a consequence of a violent dynami- 
cal event, like an interaction or a merger, followed by gas 
compression and efficient cooling. Overall, we use the cri- 
terion of considering acceptable values from 0.2 to 5 times 
the 'true' M g given by the chemical evolution code. 

3.2. An extensive grid of model template spectra 

The code allowed us to build a very large set of model 
spectra describing all possible age and mass distributions 
for the stellar populations, for the dusty ISM, and relative 
assemblies. 
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For each of the 3 histories SFR(t) we have generated 
two grids of models: one with t = 5 Myr and another with 
t = 30 Myr. Silva et al. (1998) found that the former value 
is typical for normal spirals while the latter is more suited 
for star bursting systems. Each of these grids consists of 
1400 models computed with ages to ranging from 0.2 to 
10 Gyr in steps of 0.2 Gyr and M g from to 1 (in units 
of the final mass of stars) in 28 logarithmic steps. 

In total we have therefore 1400 x 2 x 3 = 8400 
model spectra with different age, gas content, MCs es- 
cape timescale, and SFR(t) which we compared with the 
observed sample SED, allowing for the obvious scaling in 
luminosity. 

In addition we considered one further grid of spectra to 
see how our observed SEDs compare with those expected 
for spheroidal systems: for these we used the SFR(t) (c) 
model, but truncated at 3 Gyr to simulate the onset of 
a galactic wind. The adopted geometry in this case was 
a modified King profile (Eq. [3] in Silva et al. 1998) with 
r c = O.15(A//1O 11 M ) 1 / 3 Kpc. 

An example of the resulting fits to the observed broad- 
band spectra of sixsteen galaxies in our sample is reported 
in Figure |[ The analysis of the 52 fitted SEDs reveals 
the presence of two dominant different kinds of spectral 
behaviours: (a) objects which are red and show a strong 
convergence in the UV region, and (b) blue spectra that 
are flatter at all wavelengths, dominated by young stellar 
populations. 

4. RESULTS 

4-1- The space-time distribution of K-band selected 
late-type galaxies 

The grid of model spectra has been used to estimate red- 
shifts from spectral fits to the 7-band photometric data 
for the 16 galaxies lacking a spectroscopic measurement. 
As typical in cases in which such a wide spectral coverage 
(0.3 to 2.2 jum) and accurate photometry are available, 
the relative errors in z turn out to be quite small, of the 
order of ~ 10%. These broad-band spectral fits allow quite 
robust estimates of redshifts also for dusty objects, mostly 
exploiting a well-characterized feature of the optical spec- 
tra, the Balmer discontinuity, which is weakly affected by 
dust extinction (cfr discussion in section 3.2.1 in FA98). 

To check the consistency of our method, we compare 
in Figure |^ our photometric redshift predictions with the 
corresponding spectroscopic measures. The vertical error 
bars refer to different solutions at better than 90% confi- 
dence, derived from a x 2 fitting procedure using models 
(a) and (c). Fig. || shows overall good agreement within 
our modellistic uncertainties of the fits. Added to the 36 
spectroscopic redshifts, this procedure enabled us to get 
quite a reliable redshift distribution for our sample ob- 
jects. 

The distribution in redshift of a source population from 
a complete flux-limited catalogue provides a powerful con- 
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Fig. 4. Scaling of the observed average surface bright- 
ness for two faint galaxies in our sample as a function 
of redshift, according to the evolution models (a) for 
galaxy # 16 and (c) for # 50 respectively. Note that 
galaxy # 50 is the one with the lowest observed value 
of the surface brightness. The deep K-band image used 
in the primary selection has a 3er limiting brightness of 
23 K magnitudes / arcsec 2 . 



straint on its evolutionary history and formation epoch. 
This obviously assumes that we control with reasonable 
confidence all possible selection effects, in particular those 
due to the surface brightness limit, the cosmological dim- 
ming and K-corrections to the fluxes. If some morpho- 
logical criteria are at play, one needs also to understand 
how morphological appearence may evolve with redshift. 
In our case the control of the selection effects is made eas- 
ier by our primary selection in the K-band, which implies 
minimal K- and evolutionary corrections as a function of 
redshift. 

The availability of accurate measurements of the effec- 
tive radii r e allows to control the effects of the limiting 
surface brightness observable in the field in K, which we 
evaluated from the simulations described in Sect. 2 to be 
< hk >— 23 magnitudes / arcsec 2 . 

Figure ^ shows the evolution of the observed surface 
brightness for two objects in our sample (including the 
galaxy with the faintest surface brightness) as a function 
of redshift, taking into account the cosmic dimming and 
using a variety of spectral evolution patterns correspond- 
ing to the models described in Sect. 3. It is clear from the 
figure that the cutoff in surface brightness in the IRIM K 
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Min x 2 =3-9 - levels at Minx 2 + [1,5,10,15,20] Min * 2 =4 - levels at Minx 2 + [1,5,10,15,20] 
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Fig. 6. Examples of different fits to the observed broad-band spectra of object number 1, fitted with the population 
synthesis models described in Sect. 3. The figure plots \ 2 contours for the amount of dust-rich gas versus the age (top 
panels) and the corresponding best-fit spectra for two different histories of SF (lower panels). The solution on the 
right refers to the model [a] in Fig. 2 and a value for the escape timescale t a of 5 Myr (see Sect. 3.1 for details about 
these parameters). On the left: solution for model [c] and t a — 30 Myr. The contours correspond to a \ 2 increment 
of 1,5,10,15 and 20 with respect to the best fit in the grid. Therefore the second innermost contour correspond to the 
90% confidence interval. 



band image has no impact in our selection process above 
our adopted limit in total magnitude of K < 20.47, and 
the whole redshift space up to at least z — 2.5 is clearly 
accessible in principle. 

Figure [| reports the histogram of the observed red- 
shifts for our complete K < 20.47 sample (including the 
16 photometric estimates). The distribution shows pro- 
nounced peaks at z ~ 0.5 and z ~ 1, clearly indicative of 
strong inhomogeneities in the source distribution due to 
spatial clustering in the relatively small volume sampled 
by the HDFN. The uncertainty due to clustering in the 
limited volume has to be kept in mind when drawing any 
conclusions from our analysis, which require confirmation 
from surveys on more extended areas. 

A second relevant feature is apparent in Fig. ||: a cutoff 
at z ~ 1.4, with only 2 objects out of 52 found above this 
limit. 

We compared this distribution with a model prediction 
based on the local lu minosi ty function of galaxies in the 
K-band (Gardner et al. 1997 ) complemented with informa- 
tion on the contributions of various morphological classes 



from optical data (see Franceschini et al. 1998 for more de- 
tails). The luminosity function is then evolved according 
to spectral model (c), which provides a conservative esti- 
mate of the number of z > 1 galaxies (it has the minimal 
evolution rates among the three models considered). The 
other assumption we made is that the luminosity function 
changes as an effect of the evolution of the M/L ratio (M 
changes because more mass is turned into stars with time, 
L follows the evolution of the stellar populations). The 
onset of star-formation is assumed to happen at z = 4. 

As shown in the figure, these assumptions would imply 
an expected number of late- type galaxies at z > 1.5 sig- 
nificantly in excess of the observations (9 expected versus 
only 1 observed) . This result parallels a similar finding by 
FA98 and Rodighiero, Franceschini and Fasano (2000) for 
the early-type population, showing a demise of objects at 
redshifts larger than z ~ 1.3. One of the possible interpre- 
tations of this effect given by FA98, i.e. that the morpho- 
logical selection could miss galaxies with shapes deviating 
from the De Vaucouleurs profile in case of merging activity 
at these redshifts, is no longer acceptable: essentially there 
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Fig. 7. For three representative objects in our sample (# 1, 30 and 36) located at different redshift (z = 0.199, 1.355 
and 0.559 respectively) we plot various solutions for the rate of on-going SF against the average V-band effective 
extinction Ay- We report the best-fit solutions for every different history of SF considered (as described in Sect. 3.2) 
and two more corresponding models: the youngest and the oldest one with \ 2 = Xbest + <>■ The ia bels near the points 
are the \ 2 values for each solution. This figure illustrates the fact that the observed SEDs can be fitted with models 
differing in the current rate of star formation by factors up to 5-10. 



are no bright (K < 20.5) galaxies altogether at z > 1.4 in 
the HDFN area. 

4-2. Evaluating galactic ages and extinction properties 

As anticipated, if the photometric measurement of redshift 
from broad-band spectral fits is weakly affected by dust 
extinction, the estimates of most other physical parame- 
ters of gas-rich systems suffer quite more by the uncertain 
amount of dust and from the degeneracy between ages of 



stellar populations and extinction. To check this we have 
used our large grid of model spectra to study the degener- 
acy of the spectral-fitting solutions of our sample galaxies. 

Figures || plots \ 2 contours for the amount of dust- 
rich gas versus the age for two representative galaxies, as 
well as the best-fit spectra for two different histories of SF 
(model [a] and [c]). 

It is immediately apparent that, even within sets of 
models based on the same evolutionary SFR(t), a fairly 
substantial degeneracy exists between the age and amount 
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Fig. 9. Rest frame (V-K) (upper panels) and (B-J) (lower panels) colours of late-type field galaxies, compared with 
predicted values for single stellar populations with solar metallicity. The ages for the latter are indicated as well 
the mean colours of local galaxies. The left two panels refer to the observed colours, on the center we present the 
corresponding " de- reddened" colours based on our best-fit SED solutions. On the right we report the rest-frame colours 
of ellipticals and SO galaxies from FA98. Filled squares refer to objects with spectroscopic redshift, open circles to 
those with photometric redshift. The error bars shown in the central panels correspond to the uncertainty in the 
dereddening at 90% confidence. 



of dust. Furthermore, rather different star formation his- 
tories can lead to equally good fits, as detailed in Fig. || 
Figure summarizes some results of our best-fitting 
procedures for three representative objects in our sample. 
For each object, it reports various solutions for the rate 
of on-going SF and the average V-band extinction Ay, 
including the corresponding values of the x 2 ■ This figure 
illustrates the fact that the observed SEDs can be fitted 
with models differing in the current rate of star-formation 
by factors up to 5-10: a large amount of SF activity can 
be easily hidden at wavelengths below a few /im. 

Figure || details the results of two different fits to the 
observed broadband spectrum for object 30, clearly illus- 
trating the degeneracy existing between SF and extinc- 
tion. The two SEDs correspond to two solutions reported 



in Fig. 0, with values of the SFR differing by a factor ~ 2.5. 
The top panel refers to the solution 1 with x 2 = 6.77 and 
SFR— 200M Q /yr. The lower panel refers to solution 2 
with x 2 = 7.44, SFR- 400M o /yr. It is clear that if the 
analysis is confined to optical/NIR wavelenghts, it cannot 
clearly discriminate between the two solutions, whose dif- 
ferences are apparent only including the far infrared spec- 
trum, where dust re-emission would be detectable. Only 
observations of the IR spectral energy distribution, say 
between a few tenths up to a few hundreds /im, where ac- 
tively star forming galaxies emit most of the energy, would 
allow to break the present degeneracy in the solutions. 
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Fig. 5. Redshift distribution for spiral and irregular galax- 
ies brighter than K = 20.47 in the HDFN. The continuous 
line is the predicted distribution based on canonical local 
luminosity functions for late-type and irregular galaxies 
and the spectral evolution model (c) mentioned in Sect. 
3.1. 



4-3. Colours, sizes and average surface brightness of 
late-type galaxies at high redshifts 

As a first assessment of the age and extinction distribu- 
tions, we report in Figure^ the rest frame (V-K) and (B-J) 
colours as a function of redshift. As in FA98, the rest frame 
(B-J) colours are computed by interpolating the observed 
galaxy spectra using the best-fit models listed in Table 1 
(see Sect. 4.4), while the (V-K) colours require a slight 
extrapolation to longer wavelenghts. The left two panels 
refer to the observed spectra, which include the effects of 
reddening, while the colour distributions reported in the 
central panels correspond to "de-reddened" spectra (i.e. 
taking out the effect of extinction and showing the under- 
lying colour distribution) . As we see, extinction plays a sig- 
nificant role: the estimated absorption-subtracted colours 
appear on average bluer by one magnitude. De-reddened 
colours are compared with the predictions of single stellar 
populations with solar metal abundances (dashed horizon- 
tal lines). The vertical error bars in the central panels are 
the mean uncertainties related to our de-reddening proce- 
dure based on our grid of models and our adopted 90% 
confidence level. 

A comparison with the early-type galaxy sample stud- 
ied by FA98 indicates that our late-type field galax- 
ies present redder colours on average, because of extinc- 



Fig. 8. Example of two different fits to the observed 
broad-band spectrum for object number 30 in our iden- 
tification. We report the spectral galaxy emission of two 
dusty environments (Sect. 3.1): dashed line = diffuse ISM 
(cirrus), dot-dashed line = molecular clouds (MCs). The 
solid line corresponds to the total integrated spectrum of 
the galaxy. Solution 1 (top panel) and 2 (bottom panel) 
refer to that plotted in Fig. 7 for object 30, corresponding 
to a value of A v ~ 1.9, with SFR ~ 18OM /yr, X 2 = 6-77 
(for soli), and SFR - 420M Q /yr, X 2 = 7.44 (for sol2). 



tion. This evidence is stronger in the (V-K) distribution, 
where a remarkable excess of red late-types is apparent 
at z > 0.5 and V — K > 3. This illustrates that select- 
ing by colours is far more sensitive to extinction effects 
than to intrinsic differences among the stellar populations 
contributing to the flux. 

Once de-reddened, the rest-frame (B-J) colours reveal 
young stellar populations with ages from to 2 Gyrs, sig- 
nificantly bluer than those of early- type galaxies, indica- 
tive of on-going SF. The (V-K) de-reddened colours show 
a dependence on redshift: while at z < 0.6 they appear 
blue, those for galaxies at z > 0.7 are constant and quite 
red on average (2 < V — K < 3), and as red as those of 
the early- type population investigated by FA98. 

We warn that translation to age-distributions is sub- 
ject to the uncertainties in the evaluation of the effective 
extinction (see also next Section). However, the similar- 
ity in the intrinsic V-K colours of galaxies independent of 
morphology does indeed support a common age distribu- 
tion for the spheroidal stellar components in Elliptical/SOs 
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Fig. 10. Kormendy relation in the B band, i.e., the aver- 
age surface brightness versus effective radius for the bona- 
fide spirals in our high-z sample (filled cirles), compared 
with the data from a local spiral sample based on the 
RC3 (open squares). < /is > was obtained for our distant 
sources applying the K-correction and a correction for the 
internal absorption. 



and in spiral bulges, something predicted by the hierarchi- 
cal formation scenario mentioned in Section 1. 

Finally, we report in Figure [To] our measured average 
surface brightness in the B band (< /j,b >) versus ef- 
fective radius R e for the bona-fide spirals in our high-z 
sample, compared with data from a local sample based on 
the RC3. < [lb > was computed for our distant sources 
applying only the K-correction and the cosmological scal- 
ing factors. A further correction has been applied taking 
into account the effects of internal absorption. While the 
largest values of R e shown by local galaxies are missed 
by our high-z sample because we are not sampling the 
rare population of large size galaxies in the HDF limited 
space volume, there is no evidence of a significant offset 
in < jis > between the local and distant spirals within 
the large observed spreads in the data. In particular, the 
lowest surface brightness galaxies (< /is >— 22 — 23) in 
our sample are highly inclined, probably extinguished, ob- 
jects. Note that the same relation fitting the Kormendy 
relation for E/SO galaxies also fits data for spirals. 



4--4- A tentative physical characterization of late-type 
galaxies in the HDF-N 

Though aware of the uncertainties inherent in the spectral 
modelling of gas-rich systems due to the uncertain extinc- 
tion, neverthless we attempt here to estimate some basic 
physical parameters of these sources, or at least to provide 
some boundary values as found by application of our vast 
model grid. 

We report in Table 1 the formal best-fit solutions ob- 
tained from fitting the observed SEDs of our sample galax- 
ies: Ay: V band effective extinction; M to t- total baryonic 
mass divided by 10 11 solar masses; SFR: observed SFR 
in solar masses per year. 

Figure [TT|[a] plots the rate of ongoing star-formation 
SFR based on best-fit solutions versus redshift for our 
sample galaxies. The values derived in our analysis have 
a median around SFR = 10 M^fyr. Only one peculiar 
object (source number 2 in Table 1) shows an extreme 
value of SF (above ~ 500M®/yr). It is an apparently 
normal giant spiral viewed face-on, for which our spec- 
tral fit predicts a large extinction Ay = 2.2. A more 
standard extinction value {Ay = 1.3), still providing an 
acceptable fit, would still correspond to a large value of 
SFR- 2OOM /yr. 

The apparent scaling of SFR with z in Fig.|ll|[a] may 
be explained as mostly a selection effect concerning the lu- 
minousity of our objects. On the other hand, Figure [Tl|[b] 
indicates that the star formation rate is on average pro- 
portional to the intrinsic baryonic mass, such that galax- 
ies with higher SFR are typically those more massive. By 
looking at higher redshifts means to observe only the more 
luminous sources, those with larger masses. Our K-band 
selection then operates largely on the stellar mass. 

Any dependence on redshift disappears when we nor- 
malize SFR to the baryonic mass, as it is done in Figure 
|ll|[c]. The ratio M tot /SFR appearing in Fig.jljflc] gives 
the timescale for the formation of stars in our late-type 
galaxy sample. The latter does not reveal characteristics 
of violent starburst, if we consider our observed timescales 
required to convert all gas in stars: these range from 1 up 
to 20 Gyrs, and indicate a moderate star formation activ- 
ity for the present K-selected field galaxies. 

4-5. Constraints on the global star formation history: 

contributions of late-type and early-type field galaxies 

A most popular way to represent the evolutionary proper- 
ties of a population of cosmic sources is through the plot of 
the total luminosity density (or the stellar formation and 
metal production rates) in the comoving volume (Madau 
et al. 1996; Lilly et al. p.996 ). When referred to the average 
galaxy population in the field, this function was shown to 
drastically increase from the present time back to redshift 
z — 1, and to flatten off above. 
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Fig. 11. Panel a: Distribution of the on-going star forma- 
tion rate SFR versus redshift, for the best-fit solutions. 
Panel b: SFR against the total baryonic mass for best fit 
solutions. The mean uncertainty on the mass is reported. 
Panel c: Ratio of the total baryonic mass normalized to 
the on-going rate of SF for each objects (based on best-fit 
solutions). This ratio gives an estimate of the timescale 
for the conversion of gas into stars, showing a substantial 



The separate contribution of galaxies with early-type 
morphologies to the global star-formation rate per comov- 
ing volume [M© yr~ x Mpc~ 3 } has been estimated by 
FA98 using the complementary sample in the HDFN and 
population synthesis results. The outcome was that early- 
types contribute significantly to the total ^(z) mostly at 
z > 1, their fractional contribution decreasing very fast at 
lower z. 

A first reason to perform a similar computation on the 
complementary sample of spirals and irregulars is to com- 
pare the two histories of SF. A further reason of interest 
to have the full complete sample processed comes from 
recent reports claiming evidence for a more gradual de- 
cline of the galaxy ultraviolet luminosity density at z < 1 



(Cowie et al. 1999; Treyer et al. 1998), taken as an indi- 
cation of a modest evolution of the rate of SF during the 
last ~ 10 Gyrs of the galaxy cosmic history. 

An independent assessment, accounting for dust ex- 
tinction and exploiting the observed baryonic mass func- 
tion in stars through a full spectro-photometric fit to the 
SED's, would then be clearly welcome. We remember that, 
whereas this computation is relatively straightforward for 
the classified cllipticals/SO due to the lack of an ISM com- 
plicating the stellar population-synthesis fit, modelling gas 
rich late-types presents more severe problems due to the 
presence of dust. We will see later, however, that the cor- 
responding uncertainties tend to average out in the inte- 
grated form of the ^(z) function, providing a relatively 
robust result. 

We defer to the paper by Franceschini et al. (1998) 
for all details of the computation. To remind here only 
the basic steps, for all 52 objects in our complete sample 
we computed, within our grids of synthetic spectra, the 
younger and more extinguished solution. In the same way 
we determined the older solution less affected by absorp- 
tion. We computed the available comoving volumes V max 
within which the object would still be visible above the 
sample flux limit (Lilly et al. 1995, 1996). The contribu- 



tion of each galaxy to the global SFR has been estimated 
by dividing the time dependent SF rate (derived from the 
two fits) by V m ax- A correction to the comoving SF rate is 
then applied for the portion of the luminosity function not 
sampled by the present survey. Such correction is based 
on the if-band luminosity function discussed by Connolly 
et al. fll997p . The global SFR density *(t), is the summed 
contribution by all galaxies in our sample. 

The result appears in Figure [jj] in the form of the co- 
moving rate of star formation *B(z) versus redshift for the 
sample considered here (dot-dash line) , compared with the 
evolutionary path for early- type galaxies (dotted lines). 
The results in panel (a) and (b) correspond to the two ex- 
treme acceptable (at 90% of confidence) spectral solutions 
for each object, the one most extinguished and younger for 
panel (a), and the older less extinguished for panel (b). 

The disk-dominated and the irregular galaxies in the 
present sample display an evolutionary behaviour different 
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from that of bulge-dominated objects. The former appear 
to form actively stars well below z — 1, whereas the rate 
of SF for the latter is high at z > 1 but converges very fast 
at lower z. Our result for the disk and irregular galaxies 



is quite consistent with those by Brinchmann et al. (1998, 
their fig. 15), in showing a comoving SFR modestly in- 
creasing between z = 0.3 and 1. On the contrary, our re- 
sults differ significantly from Brinchmann et al. as far as 
the early- type systems are considered (in their case E/SO 
have a flat 'J'(i) in the same z-interval): we explain this as 
due to the very different procedures adopted to measure 
the function $f(z), in our case it was a global fit to the 
UV-optical-NIR SED, in their case the use of the Oil EW 
as a tracer of SF. Indeed, the latter should more likely 
trace a negligible residual of SF due to low-level merging 
activity or stellar recycling, than the global history of SF 
in these galaxies. 

Note that, despite the large uncertainties on the single 
object, the overall result is fairly well constrained between 
the two extreme solutions depicted by the shaded region 
in Figure This is due to a sort of compensation inter- 
vening in the adopted solutions: the younger-more extin- 
guished one tends to have a more intense ongoing SF activ- 
ity but less protracted in time, while the contrary happens 
for the older less-extinguished solutions. In other words, 
the baryonic mass already converted into stars and sam- 
pled by the near-IR (JHK) flux measurements as a func- 
tion of redshift, provides a more robust evaluation of the 
evolutionary SFR than the instantaneous SFR mapped 
by the short-wavelength flux. In a sense, the errorbar ap- 
pearing in Fig. [lljdoes not translates into a similarly large 
uncertainty in the prediction of Fig. [l3|, because the on- 
going rate of SF (SFR in panel [a] of Fig. [n|) and the 
timescale of SF (M to t/ SFR in panel [c]) scale inversely to 
the galaxy mass function observed at various redshifts. 

This prompted us to compare our results with those 
published by Cowie et al. (1999). This is done in Figure 
|l3| where our results appear as the shaded region, which 
is bracketed by the two solutions based on the younger- 
high extinction (upper histogram) and the older less ex- 
tinguished (lower histogram) models. The continuous line 
is a polyno mial fu nction [^(z) oc (1 + z) 1A ] quoted by 
Cowie et al. ( |l999|) as best-fitting their and Treyer's et al. 
( |l998| ) data on the time-dependent UV luminosity density. 
Within the uncertainties, our results are in quite better 
consistency with the Cowie et al. (1999) evolution ary law 
than with the dataset compiled by M adau et al. ( 1996|) , 
based on the CFRS (Lill y et a l. |l996| ) and the low-z H a 
survey by Gallego et al. (1995). 



While some discussions can be found in Cowie et al. 
( |1999| ) about possible origins for this discrepancy and on 
the consequences on this new evaluation of the evolution- 
ary SFR, we only take note here of the nice a greem ent be- 
tween our results and those of Cowie et al. ( 1999 ), based 
on quite indipendent grounds. 



It is remarkable that UV and near IR selected galaxy 
samples show such similar evolution of the comoving SFR 
density ^>(z). 

5. DISCUSSION AND CONCLUSIONS 

With the main goal to investigate systematic differences 
between early-type and late-type galaxies - as for colours, 
redshift distributions, and ages of the dominant stellar 
populations - we have analyzed a morphologically-selected 
complete sample of 52 spiral and irregular galaxies in 
the Hubble Deep Field North with total K-magnitudes 
brighter than K=20.47 and typical redshifts from z ~ 0.5 
to 1.5. The sample makes use of total photometry in the 
UBVI bands from HST and the JHK bands from ground, 
all carefully tested with an extensive set of Monte Carlo 
simulations. 

The present sample exploits in particular the ultimate 
imaging quality achieved by HST in this field, allowing 
us to disentangle among galaxy morphologies, based on 
accurate profiles of the surface brightness distributions. 

Our analysis makes also use of an exhaustive set of 
modellistic spectra accounting for a variety of physical 
and geometrical situations for the stellar populations, the 
dusty ISM, and relative assemblies. The high photomet- 
ric quality and wide spectral coverage allowed us to esti- 
mate accurate photometric redshifts for 16 objects lacking 
a spectroscopic measurement. 

We have also carefully evaluated all plausible system- 
atic effects of the selection, in particular the redshift cut- 
off implied by the limiting surface-brightness achievable in 
the reference K band image. 

A warning is in order, in any case, about the general 
conclusions derived from our sample of K-selected galax- 
ies: they should be treated with caution, due to the very 
small field of view and modest spatial sampling of the 
present survey. Ferguson et al. (2000) and Eisenhardt et 
al. (2000) estimate that the number of L* galaxies in the 
total HDF co- moving volume between z=l and z=2 is only 
a few dozens. Considering also the strong clustering in- 
ferred for Lyman break galaxies (Adelberger et al. 1998), 
statistical fluctuations imply large uncertainties on any 
conclusions based on samples like the HDF, untill more 
substantial surveys to similar depths will be made avail- 
able. 

Three the main results of our study. 

— The sample galaxies are distributed in redshift up to 
z = 1.4, but appear to be significantly missing above, 
compared with evolutionary models assuming stan- 
dard recipes for the luminosity evolution and a sub- 
stantial redshift of formation. We reported a similar 
finding in our previous study of early-type galaxies in 
the same area (Franceschini et al. 1998). Our conclu- 
sion is that, either the area has some peculiarities, or 
the underlying mass function for galaxies of all mor- 
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phological kinds has a global decline at these high red- 
shifts. Confirmation of this result will require more 
substantial sky areas to be surveyed to similar depths 
by large telescopes. 

— Differences between early- and late-types are apparent 
in the rest-frame colour distributions and the evolu- 
tionary star-formation rates per unit volume. In partic- 
ular, the short- wavelength rest-frame colours B- J, once 
dust reddening is taken into account, appear quite sig- 
nificantly bluer for late- than for early- galaxy types. 
On the contrary, the longer-wavelength colors V-K 
appear to be very similar for the two morphological 
classes. We interprete this as an indication that, while 
the stellar mix in spirals/irregulars includes young 
newly formed populations, less apparent in E/SO, the 
underlying older age component traced by the V-K col- 
ors has quite a similar origin and age distribution for 
the two galaxy categories. 

We warn, however, that the complication in spectro- 
photometric modelling introduced by dust-extinction 
in the gas-rich systems prevents to reach conclusive 
results on the source by source basis. Only future 
long-wavelength IR observations, from space (SIRTF, 
FIRST, NGST) and from ground (10 m class telescopes 
in the mid-IR and interferometers in the sub- mm), will 
allow to break down the age/extinction degeneracy. 

— We found that an integrated quantity like the 
comoving-volume star-formation rate density as a 
function of redshift $f(z) is much less affected by the 
uncertainties related to the dust distribution. The rea- 
son for this is mostly in the fact that our analysis 
is strongly constrained by the evolutionary baryonic 
mass function in stars traced by the near-IR galaxy 
luminosities, the estimate of the baryon mass at any 
redshifts being much more robust than that of the in- 
stantaneous rate of star formation (see Sect. 4.5). 

By combining this with the early-type galaxy sample 
previously studied by FA98, we find a shallower de- 
pendence of ^f(z) on z between z = 0.2 and z — 1.5 
than found by Lilly et al. (1996), i.e. *(z) cx (1 4-z) 1 ' 4 
rathen than *(z) cx (1 + z) 4 as in Lilly et al. (|1996p . 
In this redshift interval our observed ^(z) turns out 
to roughly agree with results published by Cowie et al. 
( |1999D and Treyer et al. ( |l998j ). 

Our present results, based for the first time on a care- 
ful modelling of the whole UV-optical-NIR Spectral 
Energy Distributions of galaxies, then support a revi- 
sion of the Lilly-Madau plot at low-redshifts for UV- 
and K-band selected galaxies. UV-selected and near 
infrared selected galaxy samples display a remarkably 
similar evolution of the comoving SFR density ^(z), 
at z < 1. 

The three above findings seem to favour the general 
scheme of hierarchical assembly for the formation of bright 
galaxies, envisaging their progressive build up during a 



substantial fraction of the Hubble time. After all, this is 
the most physically motivated present description. 

In this context, a warning is in order concerning some 
published specializations of the Cold Dark Matter cos- 
mogonic scheme, predicting that spheroidal galaxies in the 
field form at low redshift (z < 1) from merging of spirals 
(Kauffmann & Chariot 1995). This prediction is not sup- 
ported by our results in Fig. [l^, where Ellipticals and SOs 
appear to have been mostly formed at z = 1, whereas spi- 
rals/irregulars keep some sustained SF activity at z < 1. 
This result suggests that merging of spirals to form el- 
lipticals at low redshifts cannot be a dominant process, 
quite in ag reement with what found by Brinchmann & 
Ellis (119981) . 

It is conceivable, however, that minor modifications 
of the CDM hierachical scheme (e.g. in terms of differ- 
ent assumptions about the cosmological parameters fi m , 
Oa) can explain the observed differences between morpho- 
logical types. In our view, these differences in the pop- 
ulation histories are quite probably related to the pres- 
ence of different environments at different densities (and 
consequently different cosmic timescales of formation) in 
what we call the "field". In particular, moderately high- 
density environments (typically galaxy groups, as found 
very numerous in the spectroscopic survey by Cohen et 
al. 1999), with an accelerated cosmic timescale of evolu- 



tion and fast gas consumption, mix with truly low-density 
environments, where the transformation of primordial gas 
into stars slowly progresses during the whole Hubble time. 
We believe that the two galaxy morphologies analyzed in 
the present paper and in FA98 trace such different envi- 
ronments in the universe. 
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Fig. 12. Comoving volume star-formation rate density 
^f(t) as a function of redshift for field galaxies. The contri- 
bution of late-types to the cosmic SFR (dot dashed line) 
derived from our sample is compared in the two cases with 
the evolutionary path for early-type galaxies studied by 
FA98 (dotted line) . The solid line corresponds to the total 
amount of SF density in the field. The panels correspond 
to two different extreme solutions (see text for details) : the 
younger and more extinguished (panel [a]), the older less 
affected by dust abso rption (panel [b] ) . The data repo rted 
are from Lilly et al. ( 1996 ) and Connolly et al. ( 1997 ). 
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Table 1. Photometric data on the sample galaxies 
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18.26 


17.48 


16.79 


0.517 
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1 
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13 
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Fig. 13. The total star formation rate density per unit 
comoving volume for field galaxies (late+early-types) is 
compared with the prediction of Cowie et al. ( |1999| ), who 
found a dependence of $f(z) on z in the range 0.2 < z < 
1. of the form ^(z) oc (1 + z) 15 (continous line). The 
shaded region is bracketed by the two solutions based on 
the younger-high extinction (upper histogram) and the 
older less extinguished (lower histogram) models. See also 
caption to Fig. [l2|. 
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Abstract. We analyze a morphologically-selected com- z between z = 0.2 and z = 1.5 than found by Lilly et al 



plete sample of 52 late-type (spiral and irregular) galaxies 



in the Hubble Deep Field North with total K-magnitudes 
brighter than K=20.47 and typical redshifts z ~ 0.5 to 
1.4. This sample exploits in particular the ultimate imag- 
ing quality achieved by HST in this field, allowing us to 
clearly disentangle the early- from late-type galaxy mor- 
phologies, based on accurate profiles of the surface bright- 
ness distributions. Our purpose was to investigate system- 
atic differences between the two classes, as for colours, red- 
shift distributions and ages of the dominant stellar popu- 
lations. Our analysis makes also use of an exhaustive set 
Q of modellistic spectra accounting for a variety of physical 
and geometrical situations for the stellar populations, the 
^ dusty Interstellar Medium (ISM), and relative assemblies. 
. . The high photometric quality and wide spectral cover- 
age allow to estimate accurate photometric redshifts for 
^> 16 objects lacking a spectroscopic measurement, and al- 
low a careful evaluation of all systematics of the selection 
C« [e.g. that due to the surface-brightness limit]. This sam- 
ple appears to miss significantly galaxies above z = 1.4 
(in a similar way as an early-type galaxy sample previ- 
ously studied by us), a fact which may be explained as a 
global decline of the underlying mass function for galax- 
ies at these high redshifts. Differences between early- and 
late-types are apparent - particularly in the colour dis- 
tributions and the evolutionary star-formation (SF) rates 
per unit volume -, although the complication in spectro- 
photometric modelling introduced by dust-extinction in 
the gas-rich systems prevents us to reach conclusive re- 
sults on the single sources (only future long- wavelength IR 
observations will allow to break the age/extinction degen- 
eracy). However, we find that an integrated quantity like 
the comoving star-formation rate density as a function of 
redshift ^(z) is much less affected by these uncertainties: 
by combining this with the previously studied early-type 
galaxy sample, we find a shallower dependence of ^(z) on 



(1995). Our present results, based on a careful modelling 
of the UV-optical-NIR SED of a complete galaxy sample 
- exploiting the observed time-dependent baryonic mass 
function in stars as a constraint and attempting a first- 
order correction for dust extinction - support a revision 
of the Lilly-Madau plot at low-redshifts for both UV- and 
K-band selected samples, as suggested by independent au- 
thors (Cowie et al.|l999[). 
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1. INTRODUCTION 

Cosmogonic models, in particular those based on the hi- 
erarchical clustering of cold dark matter halos, now in- 
cluding detailed physical descriptions of gas cooling, star- 
formation and feed-back processes in the baryonic com- 
ponent, make specific predictions about the evolutionary 
history of galaxy populations as a function of their mor- 
phology. 

Basically, in the hierarchical scheme, forming galaxies 
acquire angular momentum from tidal interactions with 
the surrounding structure and then dissipate and collapse 
preferentially along the rotation vector and tend to form 
flattened rotational-supported structures (disk galaxies). 
There are indeed indications of a substantial population 
of large structures of this kind up to the highest red- 
shifts from absorption-line studies in the distant quasar 
spectra (e.g. Wolfe A.M. |1999[ ). The pressure-supported 
stellar bulges dominating E/SO galaxies, in this scheme, 
would originate from the violent relaxation and dynami- 
cal evolution following strong interactions and mergers of 
primordial disk galaxies. At the zero-th order, these mod- 
els predict that spheroidal galaxies are assembled some- 
what later than spirals, although their stellar populations 
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might not differ much if the merger occurs among gas- 
poor systems. One such extreme case has been discussed 
by Kauffmann & Chariot (1998) based on the SI = 1 stan- 
dard CDM cosmology, predicting a substantial dearth of 
spheroidal galaxies already by z=l, most of them being 
formed at lower z. 

More recently it has been pointed out that, within this 
scheme, the morphological appearence of a galaxy may 
repeatedly change with cosmic time not only from a late- 
to an early-type following a merging event, but also from a 
spheroidal to a disk configuration following the acquisition 
of new infalling gas from the environment (e.g. Ellis R.S. 



redshift. The advantage of our approach over previous at- 
tempts (Lilly et al. 1996, Cowie et al. 1999) is in our care- 



1997). This may reflect in long-wavelength (V to K) colour 
distributions un-distinguishable between early- and late- 
types, while shorter-wavelength (U to V) colours would be 
dominated by the on-going SF in disks. 

An opposite pattern is contemplated by the "tradi- 
tional" models of galaxy formation, assuming that massive 
galaxies, in particular elliptical and SO's, originated first 
at high redshifts as single entities from rapid homologous 
collapse of primordial gas. Gas-rich systems, in this view, 
form instead more quiescently from progressive inflow of 
gas into the dark matter halos during most of the Hub- 
ble time. This formation scenario then predicts a marked 
differentiation in colours and ages for the stellar popula- 
tions of the two classes of galaxies, late-type galaxies con- 
taining much younger stellar populations on average than 
the early-types. Also a substantial population of massive 
spheroids would be expected to be visible at z ~ 2 in this 
case. 

The ultra-deep integrations at various wavelengths 
performed by the Hubble Space Telescope in the Hubble 
Deep Field (Williams et al. 



1996; we consider here only 



the survey in the North area) offer an extremely valu- 
able dataset to study morphological properties of high- 
redshift galaxies. Furthermore, the very accurate photom- 
etry achievable in such deep images allows accurate es- 
timates of the photometric redshifts for vast numbers of 
faint galaxies in the field. 

We have recently exploited these data to study the 
colours, masses, age distributions, and the star-formation 
history of a sample of elliptical-SO galaxies (Franceschini 
et al. 



1998 



FA98 hereafter). The basic result was to find 
colours indicative of wide ranges of ages for the stellar pop- 
ulations and a remarkable absence of objects at z > 1.3, 
both facts telling against the predictions of the "tradi- 
tional" monolithic formation scenario. 

As a natural complement, we present in this paper an 
analysis of late-type and irregular galaxies in the HDF. 
Similarly to what we did there, our primary selection is 
in the K-band, obtained from a deep KPNO image, to 
minimize the biases in the sample due to the effects of 
K- and evolutionary corrections. The completion of our 
previous analysis of E/SO to account for the complemen- 
tary set of late-type systems is also needed for a global 
evaluation of the star-formation history as a function of 



ful treatment of dust extinction from a detailed fitting of 
the UV-optical-NIR spectral energy distributions (SED). 
In addition, the detailed knlowledge of the near-IR (NIR) 
spectrum for sources at the relevant redshifts is informa- 
tive on the baryonic mass function in stars, which provides 
an essential constraint on the cumulative star-formation 
rate as a function of time. 

In Section 2 we discuss the selection scheme and pho- 
tometric corrections used to construct a complete K-band 
flux limited sample of late-type galaxies. In Section 3 
we describe the population synthesis code that we used 
to model the optical-NIR SEDs of our sample objects, 
taking into full account the effects of a dusty interstellar 
medium in the galaxy spectra. Our main results are then 
reported in Section 4, where we perform detailed anal- 
yses of the space distributions, colours and ages of the 
stellar populations of field spiral and irregular galaxies, 
compared with ellipticals and SO. We discuss the difficul- 
ties inherent in the spectral modelling of gas-rich systems 
affected by dust extinction. We finally attempt to con- 
struct the global star formation histories of field galaxies 
(E/SO+spiral/irregulars). In Section 5 we summarize our 
main conclusions. 

We anticipate that the results of the present analysis, 
based on a survey over a very small sky field, are to be 
considered as only tentative, untill larger areas will be 
surveyed to similarly deep limits. 

We adopt Ho = 50 Km s _1 Mpc -1 throughout the 
paper. For consistency with FA98 the analysis is made 
assuming go = 0.5, and zero cosmological constant A. 

2. SAMPLE SELECTION AND PHOTOMETRY 

The Hubble Deep Field North has been observed in 4 
broad bands (F300W, F450W, F606W, F814W) for a to- 
tal of 150 HST orbits by Williams et al. (1996), and con- 
stitutes the deepest ever exposure on a small sky area. 
Accurate photometric data in the four bands have been 
published for thousands of faint galaxies by the authors. 

Dickinson et al. (1997) observed the HDF-North in 
the near-IR with the IRIM camera on the KPNO 4 m 
telescope. The camera employs a 256 x 256 NICMOS-3 
array with 0" .16 pixel -1 , but the released images were ge- 
ometrically transformed and rebinned into a 1024 x 1024 
format. IRIM exposures have been secured in the J, H 
and K filters, for a total of 12, 11.5 and 23 hours, respec- 
tively. Formal 5 a limiting magnitudes for the HDF/IRIM 
images, computed from the measured sky noise within a 
2" diameter circular aperture, are 23.45 mag at J, 22.29 
mag at H , and 21.92 mag at K , whereas the image quality 
is~l".0FWHM. 

Our sample of galaxies has been extracted from the 
HDF/IRIM A'-band image through a preliminary selec- 
tion based on the automatic photometry provided by SEx- 
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tractor (Bertin & Arnouts 1996). It is flux limited in the 
K band and it excludes early-type galaxies, i.e. objects 
whose surface brightness distribution is dominated by a 
de Vaucouleurs profile. 

To determine the limit of completeness in the K band 
(hereafter Kl) for inclusion in our sample, we followed 
the same empirical procedure described in FA98. Briefly, 
a large number of toy galaxies with exponential profiles 
(produced with the IRAF-MKOBJECTS tool) have been 
used to check the performances of Sextractor in estimating 
the K-band magnitudes (KsEx) of late-type galaxies. This 
allowed us to determine the magnitude Ku m below which 
the scatter <tk SS:c °f the measured magnitudes turns out 
to be lower than some given value <r max (the vertical line 
in Figure la corresponds to Ku m = 21 and <r max = 0.18). 
Moreover, using only galaxies with < KsEx >< Klim, we 
have derived the following empirical relation between the 
bias A A' =< KsEx > —Ktrue and the effective surface 
brightness < fsf(SEx) > (see Figure lb): 

log A A' = -1.059 + 0.188 < /if > +0.061 < > 2 

+0.006 < > 3 , (1) 

where < fif >=< fi^(SEx) > —21. This relation pro- 
vides the true total K magnitudes from the SExtrac- 
tor measured flux. Then simulations have shown that 
the K-band image has a surface brightness limit of < 
(SEx) >= 23. 

By analogy with FA98, we have first produced a cat- 
alog of morphologically selected late-type objects with 
KsEx < Ku m = 21.0; then we have used effective ra- 
dius estimates from high resolution (HST) optical imag- 
ing to derive the effective surface brightness of galaxies 
(< fif" >); finally, we have applied to the Ksex magni- 
tudes the statistical corrections given by equation [1] and 
we have included in the final sample only galaxies with 
corrected magnitudes less than or equal to Kl = Klim ~ 
AAmaj; - a max = 20.47 mag (we assume AK max ~ 0.35 
from Figure lb). 

A total of 176 objects with KsEx < 21 mag were de- 
tected by SExtractor in the IRIM A'-band image. After 
carefull inspection of the high resolution HST images, we 
rejected all elliptical and SO galaxies (including the 35 
early-type objects identified in FA98) and the stars. A few 
objects were also rejected from the sample due to their po- 
sition in the frame (at the edges of the image the noise is 
higher and the magnitude estimate is likely to be uncer- 
tain). We then produced a first preliminary, incomplete 
sample of late-type galaxies. 

The effective radii r e were estimated running SEx- 
tractor on the WFPC2 Veoe frame with the parame- 
ter FLUX -RADIUS, providing the radius containing 
half of the total emitted flux. The surface brightness 
< nf (SEx) > was evaluated for each galaxy and the sta- 
tistical corrections AA' were computed using the eq (1). 

The final complete sample of late-type galaxies with 
K C orr < 20.47 consists of 52 objects over the HDF area 
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Fig. 1. Left panel: standard deviation of the SExtractor 
magnitude estimates as a function of the average magni- 
tude for galaxies in the simulated images: the standard de- 
viation is very small for < KsEx >< 21 m .O. Right panel: 
difference between true flux and the SExtractor flux as a 
function of the average surface brightness for the subsam- 
ple of simulated galaxies with < KsEx >< 21 m .O. 



of 5.7 square arcmin. For 36 objects we have the spec 
troscopic redshift (Cohen et al 



1996, Cowie et al. 1999 



Fernandez-Soto et al. 1998| ), while for the remaining 16 we 
measured it from our photometric analysis as described 
below. 

A procedure analogous to that outlined for the K band 
was used to derive, for each object of the selected sample, 
the corrected magnitudes in the J and H bands. 

The optical magnitudes in the F300W, F450W, 
F606W and F814W bands (U, B, V and I in Table 1, re- 
spectively) have been computed again with SExtractor on 
the high resolution WFPC2 images (no corrections being 
applied in this case). Magnitudes are in the AB system, 
defined by the relation (Oke & Gunn 1983): 



AB = -2.5 log F v - 48.60 (2) 

where F v is the flux in erg cm~ 2 Hz -1 , the constant bee- 
ing choosen so that AB = V for an object with flat spec- 
trum. 

Some data on the sample are listed in Table 1. Column 
1: our identification; column 2-4: coordinates a and S (at 
J2000). To these must be added 12 hours 36 minutes (RA) 
and 62 degrees (Dec); column 5: r e is the effective radius, 
in arcsec, derived from HST images; column 6-9: optical U, 
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Fig. 2. Observed broadband spectra for sixteen galaxies in our sample, fitted with the synthetic models described in 
Section 3. 



B, V, I magnitudes in the AB system (see above); column 
10-12: near-infrared J, H, K corrected magnitudes in the 
standard system; column 13: redshift of each object. Val- 
ues in brackets are photometric redshifts, while the other 
are all spectroscopic. 



stellar medium in galaxy spectra. We defer the reader to 
that paper for a through description of this model and for 
precise definitions of the parameter, while for convenience 
we summarize the main features below. 



3. MODELLING GALAXY SEDs IN THE 
PRESENCE OF A DUSTY ISM 

The optical-NIR SEDs of our sample objects have been 
modelled u sing t he population synthesis code GRASIL 
(Silva et al. 1995 ) , taking into full account the effects (op- 
tical extinction and thermal reprocessing) of a dusty inter- 



3.1. The GRASIL code 

The code provides a self-consistent description of the for- 
mation and evolution of a galactic system in its various 
stellar and ISM components, including its secular evolu- 
tion during the Hubble time and episods of enhanced star- 
formation possibly following interactions and mergers. 
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As a preliminary step the code allows to solve the equa- 
tions ruling the chemical evolution, providing the star for- 
mation and metallicity histories SFR(t) and Z(t) as a func- 
tion of time. The computations presented here were per- 
formed adopting one-zone (no spatial dependence) open 
models including the infall of primordial gas, according to 
the standard equations of galactic chemical evolution. As 
usual, the star formation rate is determined by the amount 
of gas in the system according to a Schmidt-type law [] 

SFR{t) = vM g (t) k . 

We have generated 3 different SFR(t), in order to provide 
a wide range of spectral evolution patterns. The peak oc- 
curs at about 1, 2 and 3 Gyr (hereafter model (a), (b) 
and (c) respectively), getting broader from (a) to (c). As 
a result, half of the final stellar mass M (i.e. at 13 Gyr) 
has been assembled at galactic times of 2, 3.7 and 4.7 Gyr 
in the three cases respectively. A standard Salpeter IMF 
between 0.1 and 100 Mq is assumed. 

As described by Silva et al. (1998), GRASIL calcu- 
lates self-consistently the absorption of starlight by dust, 
the heating and thermal emission of dust grains, for an 
assumed geometrical distribution of the stars and dust, 
and a specific grain model. 

In the GRASIL model several parameters affect the 
overall modifications imprinted by dust on the SED. How- 
ever, if we confine ourselves to the attenuation of stellar ra- 
diation in the optical/UV/NIR bands, we can obtain most 
of the possible spectral behaviours by adjusting only two 
quantities: the escape timescale t (see Eq [8] in Silva et 
al. 1998 for a precise definition) of newly formed stars from 
parent molecular clouds (MCs) and the total mass of dust. 
Indeed t controls the fraction of light from very young 
stellar generations hidden inside MCs and converted to 
IR photons, since the MCs optical thickness is very high 
below ~ ljum (cfr. Silva et al. 1998). On the other hand, 
the effects of the diffuse (cirrus) dust depend on several 
quantities: the radial and vertical scale lengths for stars 
and dust distributions Rd and Zd, the residual gas in the 
galaxy M g , the dust to gas ratio S and the fraction of gas 
which is in the MCs component f mc . However we found 
that most, if not all, the possible attenuation laws of the 
diffuse dust, arising from different choices of these quan- 
tities, can be closely mimicked by simply adjusting the 
amount of gas, while fixing the other quantities to the 'typ- 
ical' values: R d = 3.5(M/1O 11 M ) 1 / 3 Kpc, z d = 0.1 R d , 
6 = 0.01 and f mc = 0.5. Obviously, while different choices 
of M g , R d , Zd, S and f mc can yield similar attenuation 
laws on the optical spectrum, the spectral shapes of the 
corresponding IR continuum re-radiation can be rather 
different . 

Strictly speaking the residual gas M g is not a param- 
eter, being instead the outcome of the chemical evolution 
code, through the Schmidt law. However we use the trick 




The code can be downloaded from http://grana.pd. astro. i 
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Fig. 3. Comparison of photometric redshifts, based on 
seven-band spectral data, with spectroscopic redshifts. Er- 
ror bars refer to interval solutions with more than 90% 
confidence level, derived from \ 2 fitting using models [a] 
and [cl. 



of forcing M g to different values, in order to describe with 
a monoparametric sequence the effects of a global attenu- 
ation on the SED. Besides this, a larger 'freedom' on M g 
takes into account that the Schmidt law should not be 
taken too literally, as a strict relationship between the to- 
tal gas content and the SFR in the system. The law may 
only provide an order of magnitude description, in par- 
ticular for the secular evolution of the SFR, the so-called 
"inactive phase" of galaxy evolution bringing essentially 
to the formation of spiral disks. Several other physical 
parameters influence the rate of star-formation with re- 
spect to the simple available amount of residual gas, in 
particular the gas pressure and temperature, which may 
drastically change as a consequence of a violent dynami- 
cal event, like an interaction or a merger, followed by gas 
compression and efficient cooling. Overall, we use the cri- 
terion of considering acceptable values from 0.2 to 5 times 
the 'true' M g given by the chemical evolution code. 

3.2. An extensive grid of model template spectra 

The code allowed us to build a very large set of model 
spectra describing all possible age and mass distributions 
for the stellar populations, for the dusty ISM, and relative 
assemblies. 
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For each of the 3 histories SFR(t) we have generated 
two grids of models: one with t = 5 Myr and another with 
t a = 30 Myr. Silva et al. (1998) found that the former value 
is typical for normal spirals while the latter is more suited 
for starbursting systems. Each of these grids consists of 
1400 models computed with ages to ranging from 0.2 to 
10 Gyr in steps of 0.2 Gyr and M g from to 1 (in units 
of the final mass of stars) in 28 logarithmic steps. 

In total we have therefore 1400 x 2 x 3 = 8400 
model spectra with different age, gas content, MCs es- 
cape timescale, and SFR(t) which we compared with the 
observed sample SED, allowing for the obvious scaling in 
luminosity. 

In addition we considered one further grid of spectra to 
see how our observed SEDs compare with those expected 
for spheroidal systems: for these we used the SFR(t) (c) 
model, but truncated at 3 Gyr to simulate the onset of 
a galactic wind. The adopted geometry in this case was 
a modified King profile (Eq. [3] in Silva et al. 1998) with 
r c = O.15(M/1O 11 M ) 1 / 3 Kpc. 

An example of the resulting fits to the observed broad- 
band spectra of sixsteen galaxies in our sample is reported 
in Figure ^[ The analysis of the 52 fitted SEDs reveals 
the presence of two dominant different kinds of spectral 
behaviours: (a) objects which are red and show a strong 
convergence in the UV region, and (b) blue spectra that 
are flatter at all wavelengths, dominated by young stellar 
populations. 

4. RESULTS 

4.1. The space-time distribution of K-band selected 
late-type galaxies 

The grid of model spectra has been used to estimate red- 
shifts from spectral fits to the 7-band photometric data 
for the 16 galaxies lacking a spectroscopic measurement. 
As typical in cases in which such a wide spectral coverage 
(0.3 to 2.2 fj,m) and accurate photometry are available, 
the relative errors in z turn out to be quite small, of the 
order of ~ 10%. These broad-band spectral fits allow quite 
robust estimates of redshifts also for dusty objects, mostly 
exploiting a well-characterized feature of the optical spec- 
tra, the Balmer discontinuity, which is weakly affected by 
dust extinction (cfr discussion in section 3.2.1 in FA98). 

To check the consistency of our method, we compare 
in Figure |3| our photometric redshift predictions with the 
corresponding spectroscopic measures. The vertical error 
bars refer to different solutions at better than 90% confi- 
dence, derived from a \ 2 fitting procedure using models 
(a) and (c). Fig. |j| shows overall good agreement within 
our modellistic uncertainties of the fits. Added to the 36 
spectroscopic redshifts, this procedure enabled us to get 
quite a reliable redshift distribution for our sample ob- 
jects. 

The distribution in redshift of a source population from 
a complete flux-limited catalogue provides a powerful con- 
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Fig. 4. Scaling of the observed average surface bright- 
ness for two faint galaxies in our sample as a function 
of redshift, according to the evolution models (a) for 
galaxy # 16 and (c) for # 50 respectively. Note that 
galaxy # 50 is the one with the lowest observed value 
of the surface brightness. The deep K-band image used 
in the primary selection has a 3<r limiting brightness of 
23 K magnitudes I arcsec 2 . 

straint on its evolutionary history and formation epoch. 
This obviously assumes that we control with reasonable 
confidence all possible selection effects, in particular those 
due to the surface brightness limit, the cosmological dim- 
ming and K-corrections to the fluxes. If some morpho- 
logical criteria are at play, one needs also to understand 
how morphological appearence may evolve with redshift. 
In our case the control of the selection effects is made eas- 
ier by our primary selection in the K-band, which implies 
minimal K- and evolutionary corrections as a function of 
redshift. 

The availability of accurate measurements of the effec- 
tive radii r e allows to control the effects of the limiting 
surface brightness observable in the field in K, which we 
evaluated from the simulations described in Sect. 2 to be 
< hk >— 23 magnitudes I arcsec 2 . 

Figure ^ shows the evolution of the observed surface 
brightness for two objects in our sample (including the 
galaxy with the faintest surface brightness) as a function 
of redshift, taking into account the cosmic dimming and 
using a variety of spectral evolution patterns correspond- 
ing to the models described in Sect. 3. It is clear from the 
figure that the cutoff in surface brightness in the IRIM K 
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Min x 2 =3-9 - levels at Minx 2 + [1,5,10,15,20] Min * 2 =4 - levels at Minx 2 + [1,5,10,15,20] 
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Fig. 6. Examples of different fits to the observed broad-band spectra of object number 1, fitted with the population 
synthesis models described in Sect. 3. The figure plots \ 2 contours for the amount of dust-rich gas versus the age (top 
panels) and the corresponding best-fit spectra for two different histories of SF (lower panels). The solution on the 
right refers to the model [a] in Fig. 2 and a value for the escape timescale t of 5 Myr (see Sect. 3.1 for details about 
these parameters). On the left: solution for model [c] and t = 30 Myr. The contours correspond to a \ 2 increment 
of 1,5,10,15 and 20 with respect to the best fit in the grid. Therefore the second innermost contour correspond to the 
90% confidence interval. 



band image has no impact in our selection process above 
our adopted limit in total magnitude of K < 20.47, and 
the whole redshift space up to at least z = 2.5 is clearly 
accessible in principle. 

Figure ^| reports the histogram of the observed red- 
shifts for our complete K < 20.47 sample (including the 
16 photometric estimates). The distribution shows pro- 
nounced peaks at z ~ 0.5 and z ~ 1, clearly indicative of 
strong inhomogeneities in the source distribution due to 
spatial clustering in the relatively small volume sampled 
by the HDFN. The uncertainty due to clustering in the 
limited volume has to be kept in mind when drawing any 
conclusions from our analysis, which require confirmation 
from surveys on more extended areas. 

A second relevant feature is apparent in Fig. [|: a cutoff 
at z ~ 1.4, with only 2 objects out of 52 found above this 
limit. 

We compared this distribution with a model prediction 
based on the local lu minosi ty function of galaxies in the 
K-band (Gardner et al. 1997 ) complemented with informa- 
tion on the contributions of various morphological classes 



from optical data (see Franceschini et al. 1998 for more de- 
tails). The luminosity function is then evolved according 
to spectral model (c), which provides a conservative esti- 
mate of the number of z > 1 galaxies (it has the minimal 
evolution rates among the three models considered). The 
other assumption we made is that the luminosity function 
changes as an effect of the evolution of the M / L ratio (M 
changes because more mass is turned into stars with time, 
L follows the evolution of the stellar populations). The 
onset of star-formation is assumed to happen at z = 4. 

As shown in the figure, these assumptions would imply 
an expected number of late-type galaxies at z > 1.5 sig- 
nificantly in excess of the observations (9 expected versus 
only 1 observed) . This result parallels a similar finding by 
FA98 and Rodighiero, Franceschini and Fasano (2000) for 
the early-type population, showing a demise of objects at 
redshifts larger than z ~ 1.3. One of the possible interpre- 
tations of this effect given by FA98, i.e. that the morpho- 
logical selection could miss galaxies with shapes deviating 
from the De Vaucouleurs profile in case of merging activity 
at these redshifts, is no longer acceptable: essentially there 
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Fig. 7. For three representative objects in our sample (^t 1, 30 and 36) located at different redshift (z = 0.199, 1.355 
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are the \ 2 values for each solution. This figure illustrates the fact that the observed SEDs can be fitted with models 
differing in the current rate of star formation by factors up to 5-10. 



are no bright (K < 20.5) galaxies altogether at z > 1.4 in 
the HDFN area. 

4.2. Evaluating galactic ages and extinction properties 

As anticipated, if the photometric measurement of redshift 
from broad-band spectral fits is weakly affected by dust 
extinction, the estimates of most other physical parame- 
ters of gas-rich systems suffer quite more by the uncertain 
amount of dust and from the degeneracy between ages of 



stellar populations and extinction. To check this we have 
used our large grid of model spectra to study the degener- 
acy of the spectral-fitting solutions of our sample galaxies. 

Figures |^ plots \ 2 contours for the amount of dust- 
rich gas versus the age for two representative galaxies, as 
well as the best-fit spectra for two different histories of SF 
(model [a] and [c]). 

It is immediately apparent that, even within sets of 
models based on the same evolutionary SFR(t), a fairly 
substantial degeneracy exists between the age and amount 
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Fig. 9. Rest frame (V-K) (upper panels) and (B-J) (lower panels) colours of late-type field galaxies, compared with 
predicted values for single stellar populations with solar metallicity. The ages for the latter are indicated as well 
the mean colours of local galaxies. The left two panels refer to the observed colours, on the center we present the 
corresponding "de-reddened" colours based on our best-fit SED solutions. On the right we report the rest-frame colours 
of ellipticals and SO galaxies from FA98. Filled squares refer to objects with spectroscopic redshift, open circles to 
those with photometric redshift. The error bars shown in the central panels correspond to the uncertainty in the 
dereddening at 90% confidence. 



of dust. Furthermore, rather different star formation his- 
tories can lead to equally good fits, as detailed in Fig. |[ 
Figure [7] summarizes some results of our best-fitting 
procedures for three representative objects in our sample. 
For each object, it reports various solutions for the rate 
of on-going SF and the average V-band extinction Ay , 
including the corresponding values of the \ 2 ■ This figure 
illustrates the fact that the observed SEDs can be fitted 
with models differing in the current rate of star-formation 
by factors up to 5-10: a large amount of SF activity can 
be easily hidden at wavelengths below a few //m. 

Figure ^| details the results of two different fits to the 
observed broadband spectrum for object 30, clearly illus- 
trating the degeneracy existing between SF and extinc- 
tion. The two SEDs correspond to two solutions reported 



in Fig. [7], with values of the SFR differing by a factor ~ 2.5. 
The top panel refers to the solution 1 with \ 2 = 6.77 and 
SFR~ 2OOM jyr. The lower panel refers to solution 2 
with x 2 = 7.44, SFR~ 4OOM /t/r. It is clear that if the 
analysis is confined to optical/NIR wavelenghts, it cannot 
clearly discriminate between the two solutions, whose dif- 
ferences are apparent only including the far infrared spec- 
trum, where dust re-emission would be detectable. Only 
observations of the IR spectral energy distribution, say 
between a few tenths up to a few hundreds //m, where ac- 
tively star forming galaxies emit most of the energy, would 
allow to break the present degeneracy in the solutions. 
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Fig. 5. Redshift distribution for spiral and irregular galax- 
ies brighter than K = 20.47 in the HDFN. The continuous 
line is the predicted distribution based on canonical local 
luminosity functions for late-type and irregular galaxies 
and the spectral evolution model (c) mentioned in Sect. 
3.1. 



4.3. Colours, sizes and average surface brightness of 
late-type galaxies at high redshifts 

As a first assessment of the age and extinction distribu- 
tions, we report in Figure I the rest frame (V-K) and (B-J) 
colours as a function of redshift. As in FA98, the rest frame 
(B-J) colours are computed by interpolating the observed 
galaxy spectra using the best-fit models listed in Table 1 
(see Sect. 4.4), while the (V-K) colours require a slight 
extrapolation to longer wavelenghts. The left two panels 
refer to the observed spectra, which include the effects of 
reddening, while the colour distributions reported in the 
central panels correspond to "de-reddened" spectra (i.e. 
taking out the effect of extinction and showing the under- 
lying colour distribution). As we see, extinction plays a sig- 
nificant role: the estimated absorption-subtracted colours 
appear on average bluer by one magnitude. De-reddened 
colours are compared with the predictions of single stellar 
populations with solar metal abundances (dashed horizon- 
tal lines). The vertical error bars in the central panels are 
the mean uncertainties related to our de-reddening proce- 
dure based on our grid of models and our adopted 90% 
confidence level. 

A comparison with the early-type galaxy sample stud- 
ied by FA98 indicates that our late-type field galax- 
ies present redder colours on average, because of extinc- 



Fig. 8. Example of two different fits to the observed 
broad-band spectrum for object number 30 in our iden- 
tification. We report the spectral galaxy emission of two 
dusty environments (Sect. 3.1): dashed line = diffuse ISM 
(cirrus), dot-dashed line = molecular clouds (MCs). The 
solid line corresponds to the total integrated spectrum of 
the galaxy. Solution 1 (top panel) and 2 (bottom panel) 
refer to that plotted in Fig. 7 for object 30, corresponding 
to a value of A v ~ 1.9, with SFR ~ 18OM /t/r, X 2 = 6.77 
(for soli), and SFR ~ 42OM /t/r, X 2 = 7.44 (for sol2). 



tion. This evidence is stronger in the (V-K) distribution, 
where a remarkable excess of red late-types is apparent 
at z > 0.5 and V — K > 3. This illustrates that select- 
ing by colours is far more sensitive to extinction effects 
than to intrinsic differences among the stellar populations 
contributing to the flux. 

Once de-reddened, the rest-frame (B-J) colours reveal 
young stellar populations with ages from to 2 Gyrs, sig- 
nificantly bluer than those of early-type galaxies, indica- 
tive of on-going SF. The (V-K) de-reddened colours show 
a dependence on redshift: while at z < 0.6 they appear 
blue, those for galaxies at z > 0.7 are constant and quite 
red on average (2 < V — K < 3), and as red as those of 
the early-type population investigated by FA98. 

We warn that translation to age-distributions is sub- 
ject to the uncertainties in the evaluation of the effective 
extinction (see also next Section). However, the similar- 
ity in the intrinsic V-K colours of galaxies independent of 
morphology does indeed support a common age distribu- 
tion for the spheroidal stellar components in Elliptical/SOs 
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Fig. 10. Kormendy relation in the B band, i.e., the aver- 
age surface brightness versus effective radius for the bona- 
fide spirals in our high-z sample (filled cirles), compared 
with the data from a local spiral sample based on the 
RC3 (open squares). < hb > was obtained for our distant 
sources applying the K-correction and a correction for the 
internal absorption. 



and in spiral bulges, something predicted by the hierarchi- 
cal formation scenario mentioned in Section 1. 

Finally, we report in Figure [lO] our measured average 
surface brightness in the B band (< hb >) versus ef- 
fective radius R e for the bona-fide spirals in our high-z 
sample, compared with data from a local sample based on 
the RC3. < hb > was computed for our distant sources 
applying only the K-correction and the cosmological scal- 
ing factors. A further correction has been applied taking 
into account the effects of internal absorption. While the 
largest values of R e shown by local galaxies are missed 
by our high-z sample because we are not sampling the 
rare population of large size galaxies in the HDF limited 
space volume, there is no evidence of a significant offset 
in < hb > between the local and distant spirals within 
the large observed spreads in the data. In particular, the 
lowest surface brightness galaxies {< Hb >— 22 — 23) in 
our sample are highly inclined, probably extinguished, ob- 
jects. Note that the same relation fitting the Kormendy 
relation for E/SO galaxies also fits data for spirals. 



4-4- A tentative physical characterization of late-type 
galaxies in the HDF-N 

Though aware of the uncertainties inherent in the spectral 
modelling of gas-rich systems due to the uncertain extinc- 
tion, neverthless we attempt here to estimate some basic 
physical parameters of these sources, or at least to provide 
some boundary values as found by application of our vast 
model grid. 

We report in Table 1 the formal best-fit solutions ob- 
tained from fitting the observed SEDs of our sample galax- 
ies: Ay: V band effective extinction; Mtot' total baryonic 
mass divided by 10 11 solar masses; SFR: observed SFR 
in solar masses per year. 

Figure |ll|[a] plots the rate of ongoing star-formation 
SFR based on best-fit solutions versus redshift for our 
sample galaxies. The values derived in our analysis have 
a median around SFR = 10 M^/yr. Only one peculiar 
object (source number 2 in Table 1) shows an extreme 
value of SF (above ~ hOOMQ/yr), It is an apparently 
normal giant spiral viewed face-on, for which our spec- 
tral fit predicts a large extinction Ay = 2.2. A more 
standard extinction value (Ay = 1.3), still providing an 
acceptable fit, would still correspond to a large value of 
SFR- 2OOM /V. 

The apparent scaling of SFR with z in Fig.pT|[a] may 
be explained as mostly a selection effect concerning the lu- 
minousity of our objects. On the other hand, Figure ^[b] 
indicates that the star formation rate is on average pro- 
portional to the intrinsic baryonic mass, such that galax- 
ies with higher SFR are typically those more massive. By 
looking at higher redshifts means to observe only the more 
luminous sources, those with larger masses. Our K-band 
selection then operates largely on the stellar mass. 

Any dependence on redshift disappears when we nor- 
malize SFR to the baryonic mass, as it is done in Figure 
fyjc]. The ratio M tot / SF R appearing in Fig.|i"l|[c] gives 
the timescale for the formation of stars in our late-type 
galaxy sample. The latter does not reveal characteristics 
of violent starburst, if we consider our observed timescales 
required to convert all gas in stars: these range from 1 up 
to 20 Gyrs, and indicate a moderate star formation activ- 
ity for the present K-selected field galaxies. 

4-5. Constraints on the global star formation history: 

contributions of late-type and early-type field galaxies 

A most popular way to represent the evolutionary proper- 
ties of a population of cosmic sources is through the plot of 
the total luminosity density (or the stellar formation and 
metal production rates) in the comoving volume (Madau 
et al. 



1996; Lilly et al. 1996). When referred to the ; 



galaxy population in the field, this function was shown to 
drastically increase from the present time back to redshift 
z ~ 1, and to flatten off above. 
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Fig. 11. Panel a: Distribution of the on-going star forma- 
tion rate SFR versus redshift, for the best-fit solutions. 
Panel b: SFR against the total baryonic mass for best fit 
solutions. The mean uncertainty on the mass is reported. 
Panel c: Ratio of the total baryonic mass normalized to 
the on-going rate of SF for each objects (based on best-fit 
solutions). This ratio gives an estimate of the timescale 
for the conversion of gas into stars, showing a substantial 



The separate contribution of galaxies with early-type 
morphologies to the global star-formation rate per comov- 
ing volume $(«) [Mq yr~ x Mpc~ 3 ] has been estimated by 
FA98 using the complementary sample in the HDFN and 
population synthesis results. The outcome was that early- 
types contribute significantly to the total ^S(z) mostly at 
z > 1, their fractional contribution decreasing very fast at 
lower z. 

A first reason to perform a similar computation on the 
complementary sample of spirals and irregulars is to com- 
pare the two histories of SF. A further reason of interest 
to have the full complete sample processed comes from 
recent reports claiming evidence for a more gradual de- 
cline of the galaxy ultraviolet luminosity density at z < 1 
(Gowie et al. 1999; Treyer et al. 1998), taken as an indi- 



cation of a modest evolution of the rate of SF during the 
last ~ 10 Gyrs of the galaxy cosmic history. 

An independent assessment, accounting for dust ex- 
tinction and exploiting the observed baryonic mass func- 
tion in stars through a full spectro-photometric fit to the 
SED's, would then be clearly welcome. We remember that, 
whereas this computation is relatively straightforward for 
the classified ellipticals/SO due to the lack of an ISM com- 
plicating the stellar population-synthesis fit, modelling gas 
rich late-types presents more severe problems due to the 
presence of dust. We will see later, however, that the cor- 
responding uncertainties tend to average out in the inte- 
grated form of the ^Sf(z) function, providing a relatively 
robust result. 

We defer to the paper by Franceschini et al. (1998) 
for all details of the computation. To remind here only 
the basic steps, for all 52 objects in our complete sample 
we computed, within our grids of synthetic spectra, the 
younger and more extinguished solution. In the same way 
we determined the older solution less affected by absorp- 
tion. We computed the available comoving volumes V max 
within which the object woul d stil l be v isible above the 
sample flux limit (Lilly et al. 1995 , 1996 ). The contribu- 
tion of each galaxy to the global SFR has been estimated 
by dividing the time dependent SF rate (derived from the 
two fits) by V max ■ A correction to the comoving SF rate is 
then applied for the portion of the luminosity function not 
sampled by the present survey. Such correction is based 
on the A'-band luminosity function discussed by Connolly 
et al. (1997). The global SFR density W(i), is the summed 
contribution by all galaxies in our sample. 

The result appears in Figure |l2] in the form of the co- 
moving rate of star formation versus redshift for the 
sample considered here (dot-dash line), compared with the 
evolutionary path for early-type galaxies (dotted lines). 
The results in panel (a) and (b) correspond to the two ex- 
treme acceptable (at 90% of confidence) spectral solutions 
for each object, the one most extinguished and younger for 
panel (a), and the older less extinguished for panel (b). 

The disk-dominated and the irregular galaxies in the 
present sample display an evolutionary behaviour different 
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from that of bulge-dominated objects. The former appear 
to form actively stars well below z = 1, whereas the rate 
of SF for the latter is high at z > 1 but converges very fast 
at lower z. Our result for the disk and irregular galaxies 
is quite consistent with those by Brinchmann et al. (1998, 
their fig. 15), in showing a comoving SFR modestly in- 
creasing between z = 0.3 and 1. On the contrary, our re- 
sults differ significantly from Brinchmann et al. as far as 
the early-type systems are considered (in their case E/SO 
have a flat ^(t) in the same z-interval): we explain this as 
due to the very different procedures adopted to measure 
the function ^(z), in our case it was a global fit to the 
UV-optical-NIR SED, in their case the use of the Oil EW 
as a tracer of SF. Indeed, the latter should more likely 
trace a negligible residual of SF due to low-level merging 
activity or stellar recycling, than the global history of SF 
in these galaxies. 

Note that, despite the large uncertainties on the single 
object, the overall result is fairly well constrained between 
the two extreme solutions depicted by the shaded region 
in Figure [D^. This is due to a sort of compensation inter- 
vening in the adopted solutions: the younger-more extin- 
guished one tends to have a more intense ongoing SF activ- 
ity but less protracted in time, while the contrary happens 
for the older less-extinguished solutions. In other words, 
the baryonic mass already converted into stars and sam- 
pled by the near-IR (JHK) flux measurements as a func- 
tion of redshift, provides a more robust evaluation of the 
evolutionary SFR than the instantaneous SFR mapped 
by the short-wavelength flux. In a sense, the errorbar ap- 
pearing in Fig. [l]does not translates into a similarly large 
uncertainty in the prediction of Fig. O, because the on- 
going rate of SF (SFR in panel [a] of Fig. |ll|) and the 
timescale of SF (M to t I SF R in panel [c]) scale inversely to 
the galaxy mass function observed at various redshifts. 

This prompted us to compare our results with those 
published by Cowie et al. (1999). This is done in Figure 
|l3| where our results appear as the shaded region, which 
is bracketed by the two solutions based on the younger- 
high extinction (upper histogram) and the older less ex- 
tinguished (lower histogram) models. The continuous line 
is a polyno mial fu nction [9(z) oc (1 + z) 1A ] quoted by 
Cowie et al. ( 1999 ) as best-fitting their and Treyer's et al. 
( 1998 ) data on the time-dependent UV luminosity density. 
Within the uncertainties, our results are in quite better 
consistency with the Cowie et al. (1999) evolutionary law 
than with the dataset compiled by Madau et al. ( 1996| ), 
based on the CFRS (Lill y et a l. |1996| ) and the low-z H« 
survey by Gallego et al. ( |1995[ ). 



While some discussions 



be found in Cc 



et al. 



(1999) about possible origins for this discrepancy and on 
the consequences on this new evaluation of the evolution- 
ary SFR, we only take note here of the nice a greem ent be- 
tween our results and those of Cowie et al. ( 1999 ), based 
on quite indipendent grounds. 



It is remarkable that UV and near IR selected galaxy 
samples show such similar evolution of the comoving SFR 
density ty(z). 

5. DISCUSSION AND CONCLUSIONS 

With the main goal to investigate systematic differences 
between early-type and late-type galaxies - as for colours, 
redshift distributions, and ages of the dominant stellar 
populations - we have analyzed a morphologically-selected 
complete sample of 52 spiral and irregular galaxies in 
the Hubble Deep Field North with total K-magnitudes 
brighter than K=20.47 and typical redshifts from z ~ 0.5 
to 1.5. The sample makes use of total photometry in the 
UBVI bands from HST and the JHK bands from ground, 
all carefully tested with an extensive set of Monte Carlo 
simulations. 

The present sample exploits in particular the ultimate 
imaging quality achieved by HST in this field, allowing 
us to disentangle among galaxy morphologies, based on 
accurate profiles of the surface brightness distributions. 

Our analysis makes also use of an exhaustive set of 
modellistic spectra accounting for a variety of physical 
and geometrical situations for the stellar populations, the 
dusty ISM, and relative assemblies. The high photomet- 
ric quality and wide spectral coverage allowed us to esti- 
mate accurate photometric redshifts for 16 objects lacking 
a spectroscopic measurement. 

We have also carefully evaluated all plausible system- 
atic effects of the selection, in particular the redshift cut- 
off implied by the limiting surface-brightness achievable in 
the reference K band image. 

A warning is in order, in any case, about the general 
conclusions derived from our sample of K-selected galax- 
ies: they should be treated with caution, due to the very 
small field of view and modest spatial sampling of the 
present survey. Ferguson et al. (2000) and Eisenhardt et 
al. (2000) estimate that the number of L* galaxies in the 
total HDF co-moving volume between z=l and z=2 is only 
a few dozens. Considering also the strong clustering in- 
ferred for Lyman break galaxies (Adelberger et al. 1998), 
statistical fluctuations imply large uncertainties on any 
conclusions based on samples like the HDF, untill more 
substantial surveys to similar depths will be made avail- 
able. 

Three the main results of our study. 

— The sample galaxies are distributed in redshift up to 
z = 1.4, but appear to be significantly missing above, 
compared with evolutionary models assuming stan- 
dard recipes for the luminosity evolution and a sub- 
stantial redshift of formation. We reported a similar 
finding in our previous study of early-type galaxies in 
the same area (Franceschini et al. 1998). Our conclu- 
sion is that, either the area has some peculiarities, or 
the underlying mass function for galaxies of all mor- 
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phological kinds has a global decline at these high red- 
shifts. Confirmation of this result will require more 
substantial sky areas to be surveyed to similar depths 
by large telescopes. 

— Differences between early- and late-types are apparent 
in the rest-frame colour distributions and the evolu- 
tionary star-formation rates per unit volume. In partic- 
ular, the short- wavelength rest-frame colours B-J, once 
dust reddening is taken into account, appear quite sig- 
nificantly bluer for late- than for early- galaxy types. 
On the contrary, the longer-wavelength colors V-K 
appear to be very similar for the two morphological 
classes. We interprete this as an indication that, while 
the stellar mix in spirals/irregulars includes young 
newly formed populations, less apparent in E/SO, the 
underlying older age component traced by the V-K col- 
ors has quite a similar origin and age distribution for 
the two galaxy categories. 

We warn, however, that the complication in spectro- 
photometric modelling introduced by dust-extinction 
in the gas-rich systems prevents to reach conclusive 
results on the source by source basis. Only future 
long-wavelength IR observations, from space (SIRTF, 
FIRST, NGST) and from ground (10 m class telescopes 
in the mid-IR and interferometers in the sub-mm) , will 
allow to break down the age/extinction degeneracy. 

— We found that an integrated quantity like the 
comoving-volume star-formation rate density as a 
function of redshift W(z) is much less affected by the 
uncertainties related to the dust distribution. The rea- 
son for this is mostly in the fact that our analysis 
is strongly constrained by the evolutionary baryonic 
mass function in stars traced by the near-IR galaxy 
luminosities, the estimate of the baryon mass at any 
redshifts being much more robust than that of the in- 
stantaneous rate of star formation (see Sect. 4.5). 

By combining this with the early-type galaxy sample 
previously studied by FA98, we find a shallower de- 
pendence of ^S(z) on z between z = 0.2 and z = 1.5 
than found by Lilly et al. (1996), i.e. *(z) oc (1 + z) 1 4 
rathen than #(z) oc (1 + z) 4 as in Lilly et al. (|1996| ). 
In this redshift interval our observed ^Sf(z) turns out 
to roughly agree with results published by Cowie et al. 
(|1999D and Treyer et al. (|1998|) . 

Our present results, based for the first time on a care- 
ful modelling of the whole UV-optical-NIR Spectral 
Energy Distributions of galaxies, then support a revi- 
sion of the Lilly-Madau plot at low-redshifts for UV- 
and K-band selected galaxies. UV-selected and near 
infrared selected galaxy samples display a remarkably 
similar evolution of the comoving SFR density ^S(z), 
at z < 1. 

The three above findings seem to favour the general 
scheme of hierarchical assembly for the formation of bright 
galaxies, envisaging their progressive build up during a 



substantial fraction of the Hubble time. After all, this is 
the most physically motivated present description. 

In this context, a warning is in order concerning some 
published specializations of the Cold Dark Matter cos- 
mogonic scheme, predicting that spheroidal galaxies in the 
field form at low redshift (z < 1) from merging of spirals 



(Kauffmann & Chariot 1998). This prediction is not sup- 
ported by our results in Fig. [l2|, where Ellipticals and SOs 
appear to have been mostly formed at z = 1, whereas spi- 
rals/irregulars keep some sustained SF activity at z < 1. 
This result suggests that merging of spirals to form el- 
lipticals at low redshifts cannot be a dominant process, 
quite in ag reement with what found by Brinchmann & 
Ellis (|l998|) . 

It is conceivable, however, that minor modifications 
of the CDM hierachical scheme (e.g. in terms of differ- 
ent assumptions about the cosmological parameters fi m , 
SI a) can explain the observed differences between morpho- 
logical types. In our view, these differences in the pop- 
ulation histories are quite probably related to the pres- 
ence of different environments at different densities (and 
consequently different cosmic timescales of formation) in 
what we call the "field". In particular, moderately high- 
density environments (typically galaxy groups, as found 
very numerous in the spectroscopic survey by Cohen et 
al. 1999), with an accelerated cosmic timescale of evolu- 



tion and fast gas consumption, mix with truly low-density 
environments, where the transformation of primordial gas 
into stars slowly progresses during the whole Hubble time. 
We believe that the two galaxy morphologies analyzed in 
the present paper and in FA98 trace such different envi- 
ronments in the universe. 
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Fig. 12. Comoving volume star-formation rate density 
^(t) as a function of redshift for field galaxies. The contri- 
bution of late-types to the cosmic SFR (dot dashed line) 
derived from our sample is compared in the two cases with 
the evolutionary path for early-type galaxies studied by 
FA98 (dotted line). The solid line corresponds to the total 
amount of SF density in the field. The panels correspond 
to two different extreme solutions (see text for details): the 
younger and more extinguished (panel [a]), the older less 
affected by dust abso rption (panel [b]). The data repo rted 
are from Lilly et al. ( 1996 ) and Connolly et al. ( 1997 ). 
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Table 1. Photometric data on the sample galaxies 
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Fig. 13. The total star formation rate density per unit 
comoving volume for field galaxies (late+early-types) is 
compared with the prediction of Cowie et al. ( 1999| ), who 



found a dependence of ^(z) on z in the range 0.2 < z < 
1. of the form $(z) oc (1 + z) 1 5 (continous line). The 
shaded region is bracketed by the two solutions based on 
the younger-high extinction (upper histogram) and the 
older less extinguished (lower histogram) models. See also 
caption to Fig. M3. 



